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Abstract. This article investigates the computation of the eigenmodes of 
the Laplacian operator in multi-connected three-dimensional spherical spaces. 
General mathematical results and analytical solutions for lens and prism spaces 
are presented. Three complementary numerical methods are developed and 
compared with our analytic results and previous investigations. The cosmological 
applications of these results are discussed, focusing on the cosmic microwave 
background (CMB) anisotropics. In particular, whereas in the Euclidean case 
too small universes are excluded by present CMB data, in the spherical case there 
will always exist candidate topologies even if the total energy density parameter 
of the universe is very close to unity. 



PACS numbers: 98.80.-q, 04.20.-q, 02.040.Pc 

1. Introduction 

The search for the topology of our universe has made tremendous progress in the past 
years and several methods have been designed using either galaxy catalogs or the CMB 
(see e.g. [[[], || for general reviews). The most promising dataset that can eventually 
contain a topological signature is the cosmic microwave background (CMB) in the 
form of pattern correlations (such as homologous circles in the sky ||] , or anomalously 
large temperature correlations in a set of directions Q , see || for a recent review on 
the CMB methods) or non Gaussianity ||. 

The detectability of the topology in datasets such as those that will be made 
available by the MAP and Planck || satellite missions requires to simulate maps 
with the topological signature for a large set of topologies. These maps will have 
mainly two uses: first, they will allow us to test the detection method and for instance 
estimate its running time and second, once all sources of noises are added, it will help 
us investigating to which extent a given method is well suited to detect the topological 
signal and indeed if it is not blurred (in the same spirit as the investigation of the 
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"crystallographic" methods based on galaxy catalogs Q). A prerequisite for any 
further study is thus to simulate CMB maps with a topological signal. 

At present, the status of the constraint on the topology of the universe is sparse. 
Concerning locally Euclidean spaces, it was shown on the basis of the COBE data that 
the size of the fundamental domain of a 3-torus has to be larger than L > 4800 ft." 1 Mpc 
[ pT| |ll|, [l^, |l3| . This constraint does not exclude a toroidal universe since there can 
be up to N — 8 copies of the fundamental cell within our horizon. This constraint 
relies mainly on the fact that the smallest wavenumber is 2tt/L, which induces a 
suppression of fluctuations on scales beyond the size of the fundamental domain. This 
result holds only for the case of a vanishing cosmological constant and was generalized 
to all Euclidean manifolds JlJ] . A non- vanishing cosmological constant induces larger 
scale cosmological perturbations, via the integrated Sachs- Wolfe effect. For instance 
if Oa = 0.9 and f2 m = 0.1, the former is relaxed to allow for N = 49 copies of the 
fundamental cell within our horizon. This constraint is also milder in the case of 
compact hyperbolic manifolds and it was shown |l7j that the angular power 

spectrum was consistent with the COBE data on multipoles ranging from 2 to 20 for 
the Weeks and Thurston manifolds. Another approach was developed in |l8|, [H], 
and is based on the method of images. Only one spherical space using this method of 
images was considered in the literature, namely the case of the projective space 

Note that multiconnectedness breaks global homogeneity and isotropy (except 
for the particular case of the projective space). It follows that the temperature 
angular correlation function C will depend on the position of the observer and on 
the orientation of the manifold, which is at odd with the standard lore. In a simply 
connected universe the angular correlation function depends only on the angle between 
the two directions whereas in a multi-connected universe, it will depend on the 
two directions. It follows that the coefficients Ci of the decomposition of C into 
Legendre polynomials, obtained by averaging over the sky, loose much of topological 
information. As clearly explained in ||, C can be decomposed into an isotropic and 
an anisotropic part and the Ci depend solely on the former. The Ci alone are a poor 
indicator of the topology, despite the fact that they can help constraining the topology, 
which backs up the necessity to study the full sky map. 

In standard relativistic cosmology, the universe is described by a Friedmann- 
Lemaitre spacetime with locally isotropic and homogeneous spatial sections. These 
spatial sections can be defined as the constant density or time hypersurfaces. With 
such a splitting, the equations of evolution of the matter and geometry perturbations 
that will give birth to the large scale structures of the universe reduce to a set of 
coupled differential equations involving a Laplacian (see e.g. Q). This system is 
conveniently (numerically) solved in Fourier space but this requires to determine 
the eigenmodes and eigenvalues of the Laplacian through the generalized Helmoltz 
equation 

AV q = W (1) 

The Laplacian in Eq. (|l|) is defined as A = DiD 1 , Di being the covariant derivative 
associated with the metric 7^ of the spatial sections («, j — 1..3). The eigenmodes on 
which any function can be developed encode the boundary conditions imposed by the 
topology, and any function developed on this basis will satisfy the required boundary 
conditions. 

Concerning the computation of the eigenmodes, the case of Euclidean manifolds 
can be solved analytically and many numerical investigations of compact hyperbolic 
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manifolds have been performed |l6|, |24|, ^5|, 26 . The case of spherical manifolds 
has been disgarded during a long time. The recent measurements of the density 
parameters let the room for our universe to be slightly positively curved since they 
are estimated to lie in the range Slo = ^m + ^A = I- 1 1-0.12 t° 95% confidence p7| . 

The goal of this article is to focus on the computation of the eigenmodes of the 
Laplacian in spherical spaces, having in mind their use to simulate CMB maps. In 
Section & we recall the main mathematical results on the classification of spherical 
spaces p8j and some analytical results on the eigenmodes of the Laplacian such as 
the determination of the smallest wavenumber and its multiplicity. In section [jj we 
describe the three numerical methods that are used to compute the eigenvalues and 
eigenmodes and we discuss their precision and efficiency. We then describe briefly, 
in Section f|, a method to obtain these eigenfunctions analytically in the particular 
cases of lens and prism spaces. Such a result is of importance to compare with the 
output of numerical computations. After describing some statistical properties of the 
eigenmodes in Section |^, we present, in Section ^, some cosmological consequences 
of our result. We estimate the effect of the topology on the Sachs- Wolfe plateau. 



Appendix A gathers most of the results on the eigenfunctions of the Laplacian in 



homogeneous and isotropic 3-dimensional simply connected spaces. 

Notations: The local geometry of the universe is described by a Friedmann-Lemaitre 
metric 

d S 2 = -dt 2 + a 2 (i)[dx 2 + / 2 (x)dfi 2 ] (2) 

with f(x) = (sinhx, x, sm x) respectively for hyperbolic, Euclidean and spherical 
spatial sections. a(t) is the scale factor, t the cosmic time and dfi 2 = d9 2 + sin 2 9dip 2 
the infinitesimal solid angle, x is t ne (dimensionless) comoving radial distance in units 
of the curvature radius Rc ■ 

We use the embedding of the 3-sphere S 3 in 4-dimensional Euclidean space 
by introducing the set of coordinates (x M )n=o..3 related to the intrinsic comoving 
coordinates (x,#, <p) through (see e.g. Ref. pH) 

x = cos x 

xi = sin x sin 9 sin tp 

X2 = sin x sin 9 cos tp 

X3 = sin x cos 9, (3) 

with 0<x< 7r I 0<#<7r and < tp < 2ir. The 3-sphere is then the submanifold of 
equation 

x^Xfj. = xl + x\ +x 2 2 +xj = +1, (4) 

where x M = S^x 1, '. The comoving spatial distance d between any two points x and y 
on S 3 can be computed using the inner product x^y^ and is given by 

d[x, y] = arc cos [x^y^] . (5) 

The volume enclosed by a sphere of radius x is, in units of the curvature radius, 

Vol(x) =^(2 X -sin2 X ). (6) 

We consider 3-dimensional multi-connected manifolds of the form — S 3 /F 
where the holonomy group|] F is a discrete subgroup of 5*0(4) that acts without fixed 

t We use this term to mean the geometric version of the group of covering transformations, that is, 
elements of the holonomy group are isometries, while covering transformations are homeomorphisms, 
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point on the 3-dimensional covering space S 3 . It is isomorphic to the first fondamental 
group ir 1 (A43). The elements g £ T are isometries and can be expressed as 4 x 4 
matrices acting on the elements of the 4-dimensional embedding Euclidean space. 

The transformation (|^) between the intrinsic coordinates (x, 0,4>) and x^ is 
bijective. In the following of the article, if / is a function on 5 3 , we will loosely 
write f{x) for /(x, 8, (f) and for instance f(gx) will refer to /(%', 0' , <f)') with (x' , 0' , 0') 
being the intrinsic coordinates of the image x' = gx obtained from Eq. (||) [see e.g. 
Eq. (A. 14) for an example]. 

To finish, let us remark that any eigenmode of a multi-connected manifold S 3 /T 
lifts to a r-invariant|] eigenmode of S 3 , and conversely each T-invariant eigenmode of 
S 3 projects down to an eigenmode of S 3 /T. Thus, with a slight abuse of terminology, 
we may say that the eigenmodes of S 3 /T are the T-invariant eigenmodes of S 3 . 



2. Mathematical results on the Laplacian operator in spherical spaces 

2.1. Classification of spherical manifolds 

In our preceding article ]2S| ], we presented in a pedestrian way the complete 
classification of 3-dimensional spherical topologies and we described how to compute 
their holonomy transformations. We briefly outline the main results of this study, 
mainly to set our notations. The classification of constant curvature spherical 
manifolds was first presented in . 

The isometry group of the 3-sphere is G = 50(4) and one needs to determine 
all the finite subgroups of 50(4) acting fixed point freely. Every isometry of 50(4) 
can be uniquely decomposed as the product of a right-handed (R) and a left-handed 
(L) Clifford translation!], up to a multiplication by -1 of both factors. Besides, S 3 
enjoys a group structure as the set S 3 of unit length quaternions. It can then be 
shown that each right-handed (resp. left-handed) Clifford translation corresponds to 
a left (resp. right) multiplication of S 3 [q — > xq (resp. q — > qx)] so that the two 
groups of right-handed and left-handed Clifford translations are isomorphic to S 3 . It 
follows that 50(4) is isomorphic to S 3 x 5 3 /{±(l, 1)} so that the classification of the 
subgroups of 50(4) can be deduced from the classification of all subgroups of S 3 . It 
can then be shown that there exists a two-to-one homomorphism from S 3 to 50(3) 
the subgroups of which are known. It follows that the finite subgroups of S 3 are: 

• the cyclic groups Z n of order n, 

• the binary dihedral groups D* m of order 4m, m > 2, 

• the binary tetrahedral group T* of order 24, 

• the binary octahedral group O* of order 48, 

• the binary icosahedral group I* of order 120, 

where a binary group is the two-fold cover of the corresponding group. 

^From this classification, it can be shown that there are three categories of 
spherical 3-manifolds. 

• The single action manifolds are those for which a subgroup R of S 3 acts as pure 
right-handed Clifford translations. They are thus the simplest spherical manifolds 
and can all be written as S 3 /T with T = Z n , D* m , T*,0*,I*. 

§ A T-invariant function / is a function satisfying f(x) = f(gx) for all g 6 T and for all x. 
| Clifford translations are isometries that translate all points the same distance. 
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• The double-action manifolds are those for which subgroups R and L of S 3 act 
simultaneously as right- and left-handed Clifford translations, and every element 
of R occurs with every element of L. There are obtained for the groups r = Ti x T 2 
with (ri,r 2 ) = (Z m , Z n ), (D* m , Z n ), (T*,Z n ), (0*,Z n ), (I*,Z n ) respectively with 
gcd(m,n) = 1, gcd(4m,n) = 1, gcd(24, n) = 1, gcd(48, n) = 1, gcd(120, n) = 1. 

• The linked-action manifolds are similar to the double action manifolds, except 
that each element of R occurs with only some of the elements of L. 

The classification of these manifolds is summarized in the figure 8 of . 

We also define a lens space L(p, q) by identifying the lower surface of a lens-shaped 
solid to the upper surface with a 2nq/p rotation for relatively prime integers p and q 
with < q < p. Furthermore, we may restrict our attention to < q < p/2 because for 
values of q in the range p/2 < q < p the twist 2nq/p is the same as — 2tt(p — q)/p, thus 
L(p, q) is the mirror image of L(p,p — q). Lens spaces can be of any of the category 
described above and their classification is detailed in figure 9 of 

2. 2. Spectrum of spherical manifolds 

The Helmoltz equation (|l|) is usually rewritten in terms of the quantity (3 defined as 

P 2 = Q 2 + 1 (7) 

and takes the form 

A¥ M =-k(fc + 2)¥ M , (8) 

after using the change of variable /3 = k + 1 and where s is an integer labelling the 
modes of same eigenvalue (see Ref. j32| for the properties of the Laplacian operator). 
It follows that the eigenvalues of the Laplacian on S 3 are k(k + 2), k being an integer, 



and that the multiplicity of each eigenvalue is (k + l) 2 [see Appendix A for details]. 
The two modes k — and k = 1 are gauge modes since they respectively correspond to 
a change in the curvature radius (homogeneous deformation) and in a displacement of 
the center of the 3-sphere and are thus physically not relevant |33| . It is thus clear that 
the eigenvalues of the Laplacian on a multi-connected spherical manifold = S 3 /F 
will be a subset of the eigenvalues k(k + 2) of the Laplacian on S 3 . 

For any spherical manifold AAz, we can introduce its (discrete) spectrum as the 
set of all eigenvalues of the Laplacian 

S P (X 3 ) = {0 = A < A x < A 2 < . . . < Ai < . . .} (9) 

and two Riemannian manifolds 7W3 and M'3 are said to be isospectral if Sp(A^3) = 
SpCM' 3 ). 

Ikeda and Yamamoto J34j studied the spectra of 3-dimensional lens spaces and 
demonstrated that if two 3-dimcnsional lens spaces with fundamental group of order q 
are isospectral to each other, then they are isometric to each other. Ikeda |3f| extended 
this result to show that if two 3-dimensional spherical manifolds are isospectral 
then they are isometric. Combining this with previous results, it follows that a 3- 
dimensional spherical manifold is completely determined, as a Riemannian manifold, 
by its spectrum. Higher-dimensional generalizations appear in |3^| . 

More related to our purpose is the work by Ikeda J37| in which the spectra of 
single-action manifolds are determined. This result is of great interest for comparison 
with our numerical computation. Unfortunately, it will give us only the wavenumbers 
with their multiplicity for a restricted set of topologies and it does not determine the 
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Manifold 



Eigenvalue k(k + 2) 



multiplicity 



S 3 /Z2p 
S 3 /Z2p+1 

S 3 /D* m 
S 3 /T* 

s 3 /o* 
s 3 /r 



k > even 

k > 2p + 1 odd 

k > even 

k even 

k > m odd 

^ 1,2,5 
fc ^ 1,2,3,5,7,11 
k ^ 1,2,3,4,5,7,8, 

9,11,13,14,17,19,23,29 



( k + !)Ei=o n fe' 

+ 1 )Ef=o n ^ 

(2fc+ l)([fc/m] + 1) 

(2fc+ l)[fc/m] _ 

(2k + 1)(1 + 2[fc/3] + [fc/2] - fc) 

(2fc + 1)(1 + [fe/4] + [fc/3] + [fc/2] - fc) 

(2fc + 1)(1 + [Is/ 5] + [fc/3] + [fc/2] - fc) 



Table 1. Spectra of single-action manifolds. For spaces other than S 3 /Z n we 
have set k = 2k with k being an integer, [p] refers to the integer part of p andnki 
is defined by n^j = 1 if k = 2£(modra) and otherwise. Adapted from Ref. [B7J. 



Manifold 

S 3 /Z n , n>2 
S 3 /D* 2 

S 3 /D; n , m > 2 
S 3 /T* 

s 3 /o* 
s 3 /r 



First eigenvalue fc(fc + 2) 



8 
8 

24 
24 
48 
80 
168 



multiplicity 



2 


9 


2 


3 


4 


10 


4 


5 


6 


7 


8 


9 


12 


13 



Table 2. Value and multiplicity of the first non zero eigenvalue for the single- 
action manifolds. From Ref. |37| . Note that we have corrected a mistake on the 
multiplicity of the first eigenvalue of S 3 / D*. 



eigenfunctions. In tables |l| and [2| we sum up the main results on the wavenumbers 
of single-action manifolds. In figure |l|, we compare the spectra of the projective space 
P 3 = S 3 jZi and of the 3-sphere. Due to the Z^ symmetry, half of the modes are 
lost because they are not invariant under the antipodal map. In figure || we compare 
the spectra of the single- action manifolds obtained from the groups Z7, Z$, Z g , 
-D3, and D\ and figure || compares the ones obtained from the binary tetrahedral, 
octahedral and icosahedral groups. 

3. Numerical determination of the eigenmodes 

To determine the eigenvalues and eigenmodes of the Laplacian, one has to find a way 
to take into account the boundary conditions imposed by the topology. Different 
routes have been investigated. The case of locally Euclidean manifolds is somehow 
trivial since the problem can be solved analytically (see e.g. p8|]). The hyperbolic case 
was first addressed using the boundary element method first developed by Aurich and 
Steiner |^3| for the study of 2-dimensional hyperbolic surfaces. Inoue |23] developed 
the direct boundary element method and was the first to determine precise eigenmodes 
of 3-dimensional compact hyperbolic manifolds and get the 36 first eigenmodes of 
Thurston space for fc < 10 0] and then for fc < 13 [[16) (see also Ref. for the 

first computation of the eigenmodes of a cusp manifold). Recently a new method was 
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eigenvalue k(k+2) 



Figure 1. Comparison of the spectra of the 3-sphere and of the projective space. 
Half of the modes are lost due to the reflection symmetry (since antipodal points 
are identified, the eigenfunction must vanish on the equator). 
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Figure 2. Spectra of single-action manifolds generated from [left] the cyclic 
groups Z7, Zs and Zg and from [right] the binary dihedral groups D^-, D$ and 
D* 



proposed by Cornish and Spergel |?6j in the framework of hyperbolic spaces; such a 
method can be adapted to the case of spherical spaces, and thus will be described in 
more details below as the "ghosts method" . 

In this section, we present three independent methods to compute the eigenmodes. 

3.1. Ghosts method 

The ghosts method is based on the idea that any square integrable function in 
L 2 (X/T), X being the universal covering space, satisfies 

tf(af) = V(gx) (10) 

for all g £ T. Any function of L 2 (X/T) can be lifted to a (T-invariant) function of 
L 2 (X) and, reciprocally, any T-invariant function projects down to a function on X/T. 
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Figure 3. Spectra of single-action manifolds generated from the binary 
tetrahedral, octahedral and icosahedral groups. 



It follows that the eigenmodes of the Laplacian can be decomposed as 

^0,a( x )=^2 ZpMmypimtx) (H) 

e=o m=-e 

where there is no summation on /? because the eigenfunctions of the universal 
covering space, y^ m = Rf3i(x)Y£ m (0,4>) ( see Appendix A ) form a complete and 
linearly independent family. Note that the coefficients £p tS em are obtained once the 
holonomies are known and they will depend on the base point. 

Now, choose randomly d points, Xi, in the fundamental domain and consider the 
m images of each point up to a distance p max . We also need to truncate the sum ( pi] ) 
to a maximum value L of I. Each point generates n^n,; + l)/2 constraints of the form 



(12) 



^(i.A9aXi) - Wp^ljbXi) = 

for a ^ b and a, b = . . . \T\. With the decomposition ( pi] ) the set of constraints (|12| 
for all the d points Xi takes the form 



A-pVp, s = 

where the (L + l) 2 -components vector vp tS is defined by 



Vj3a 



£,P,sLL 



(13) 



(14) 



The matrix Ap with M = J2j=i n j( n j + -0/2 rows and N = (L + l) 2 columns i 
defined by^j| 

yp00(j9lXl) ~ 3^00(32^1) ■ • • ypLL{,9lXl) - y/3LL(92Xl) 
% Note that there is a typo in the equation (2.3) of |26|. 



yf3LL(g a Xd) - ypLL(gbXd) 



(15) 
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If M > N the system (|TJ) is over constrained and has a solution if and only if q 
is an eigenvalue of the compact space. We named this method as the ghosts method 
because of its close analogy with the simulation of galaxy catalogs in a multi-connected 
universe @ || g 

Numerically, one uses a single value decomposition (SVD) method to extract the 
eigenmodes, i.e. which form a basis of Ker(A l g). This decomposition is based on the 
theorem stating that any M x N matrix A with M > N can be written as the product 

A = UD*V (16) 

where U is a M x N unitary matrix, D = diag(w l7 . . . ,u>n) is a N x N diagonal 
matrix, the Wi being non negative and V is a N x ./V unitary matrix. The columns 
of U associated to non-zero if, form an orthonormal basis spanning the range of A^. 
The columns of V corresponding to the vanishing wi form a basis of the nullspace of 
A^ since it is solution of the equation (fl3|), i.e. of the subspace Ep ] of the eigenmodes 
of X/T with eigenvalue q. 

This method was first implemented |2^] to compute the lowest eigenvalues and 
eigenfunctions for 12 hyperbolic manifolds. We have extended the method to include 
any topology. In hyperbolic spaces, we recover the results of ]2q| and in the Euclidean 
case we compared our result with the analytic solutions. When the universal covering 
space is not compact, one has to specify the value of the degree of constraint c = M/N, 
of the cut-off L and of the maximum radius p max up to which the images are considered. 
Our main interest in this article is the case of spherical spaces and it enjoys a number 
of simplifications. First, we can set p max = tt since the group T is finite. It follows 
that any of the d points has exactly |T| — 1 images so that M — d|r|(|r| — l)/2 and 
we end with only two free parameters (L, d) or equivalently (L,c). This simplifies the 
discussion concerning the choice of the different cut-offs L. 

The parameter c is adjusted numerically. It needs to increase fast with the order 
|r| of the group, but for each order c needs to increase slowly with the increase of (3. 
For example for |T| = 8, c = |T| + j3 works very well until (3 = 15. For (3 > 15, we 
need either a smaller c or a c that increases slower with (3. The constraint determines 
the d points that we need to generate randomly in the fundamental domain for each 
(3. Apparently the only reason for choosing a specific value for d is that the method 
does not work without a good balance between M (number of rows) and N (number 
of columns) for each |T| and (3. A wrong choice of c could cause the failure of the 
method. The results of this method have been compared with Ikeda's results and 
agree with them. 



3.2. Averaging method 

When focusing on spherical spaces, one can take into account the fact that the number 
of images of any point is exactly given by the order of the group, and is thus finite, 
to develop another numerical method for computing the eigenfunctions, completely 
independent from the previous one. 

This method is based on the two remarks that 

• if ^ is an eigenmode, then ^ o g is also an eigenmode 

• for any function / on S 3 , then / defined by 
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is a T-invariant function of L^(S 3 ), the indice recalling the T-invariance, which 
can thus be identified to a function of L 2 (S 3 /T). The operation 

avr: f ±f (18) 

is a projection on the subspace of T-invariant functions (since av 2 , = 1 and 
av r(/) = / iff / is T-invariant). This induces an equivalence relation g ~ / 
iff avr(.g) = avr(/) so that L 2 (S 3 /T) is just given as the set of the functions of 
L 2 (S 3 ) modulo the equivalence relation, i.e. {avr(/), / G L 2 (S 3 )} 

It follows that the eigenmodes of the Laplacian on S 3 /T are explicitly given in terms 
of the eigenmodes of the Laplacian on the universal covering space (see Appendix A ) 

by 

1 1 ger 

Indeed, the (k + l) 2 linearly independent eigenmodes yum project down to [k + l) 2 

[r] 

T-invariant eigenmodes W fe . This set of functions is a generator of the eigenspace 
e\ 1 but it is not a free family because we expect dim(i^ ')<(&+ l) 2 . We will thus 
have to pick up the dim(E [ k ] ) independent functions of the family (|l9|). This can 
be performed by using an orthonormalisation procedure (the classical Gram-Schmidt 
method itself being numerically disastrous). 

Numerically, we compute the average family ( |l9| ) and then decompose it on the 
basis y k em as 



^ [ ] .s — £p, s e m y pem (20) 



and use the SVD method (as described in the previous section) to perform the 
orthonormalisation. For that purpose we have to choose d = (k + l) 2 and L = (3 
so that there is no free parameter to choose and 

e -5(sT5)'™n-i)- (^1) 



3.3. Projection method 

The ghosts method and the averaging method, presented in the two previous sections, 
both compute an orthonormal basis for the space of eigenmodes of S 3 /T, so that 

we may later compute a random eigenmode 'I'^ 1 of S" 3 /r as a linear combination (|ll|). 

[r] 

The projection method, by contrast, computes the random eigenmode ^ ^ directly, 
without explicitly computing a basis for the eigenspace. Throughout this section 
we assume the wavenumber k is fixed, thus fixing the eigenvalue k(k + 2) and the 
parameter (3 — k + 1 as well. 



3.3.1. The algorithm The idea is to first compute a random eigenmode of S 3 , and 
then project that random eigenmode down to an eigenmode of S 3 /T. More precisely, 
one proceeds as follows: 
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Step 1: Construct a random eigenmode on S 3 . 

Construct a random eigenmode on S 3 as a linear combination 

*/3 = ^ (t3£ m ypim- (22) 
l,m 

Choose the coefficients Cp£m relative to a Gaussian distribution with mean and 
standard deviation 1. The resulting point ((p£ m ) G R* 3 will, in effect, be chosen 
relative to a spherically symmetric distribution, because the product measure 

nexp(-Cj, m /2)=exp(-r 2 /2) 

depends only on the radial distance r 2 = Y^z m (pem m ■ Note that the 
expected value E((pc m ) of each coefficient (pe m is 1, so the expected value of the 
squared radius r 2 is simply the dimension (3 2 of the space: 

E(r 2 ) = E(£ C} lm ) = J2 E(Chm) = £ 1 = 

In the infinitely unlikely event that r = 0, discard the randomly chosen Cpe m and 
choose a new set. 

Step 2. Construct the average &vr(^p). 

Given the eigenmode <J p of S 3 , define the T-invariant eigenmode avr(^p) of S 3 /Y 
via the averaging formula ([l7]). As mentioned earlier, a T-invariant eigenmode of 
S 3 corresponds to an eigenmode of the quotient space S 3 /Y. In principle we can 
evaluate formula (jf^) for any igS 3 , but in practice we need evaluate it only for 
those x lying on the last scattering surface. 

Computational note: We begin with x in rectangular coordinates (xi, X2, £3, £4), 
so that g{x) may be computed quickly and easily as a 4 x 4 matrix times a 4- 
clcment vector. We then convert the result to spherical or toroidal coordinates 
(x, 0, <p) for efficient evaluation of Eq. (p2|). 

3.3.2. Proof of orthogonality Let n — f3 2 be the dimension of the full /3-eigenspace 
of S* 3 , and let m be the (unknown) dimension of the T- invariant subspace, i.e. m 
is the dimension of the /3-eigenspace of S 3 /Y. Thus the averaging operator avr() in 
Step 2 projects the eigenspace R n of S 3 down onto the subspace R m of T-invariant 
eigenmodes. The question is, relative to the standard inner product in function space, 
is this an orthogonal projection from R" to R m ? If so, then a spherically symmetric 
distribution of points in R" (corresponding to random eigenmodes of S 3 ) will project 
to a spherically symmetric distribution of points in R m (corresponding to random 
eigenmodes of 5 3 /r). However, if the projection is not orthogonal - for example if 
it involves a shearing motion - then a spherically symmetric distribution of points 
in R" will not project to a spherically symmetric distribution in R m , but rather to 
a distribution that is slanted to one side or another. Fortunately the projection is 
orthogonal, as shown in Proposition I. 

Proposition I. Let H be the space of all (3-eigenmodes on S 3 , and Hp be the space of 
(5-eigenmodes that are invariant under the action of a finite group Y C 50(4). Then 
the averaging function avr defined in Eq. mA) maps H orthogonally onto Hr ■ 
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Proof. Let H± be the orthogonal complement of Hy in H. That is, let H± consist of 
those eigenmodes in H that are orthogonal to all eigenmodes in Hp. Thus H ~ R™, 
Hp ~ R m , and H± :_ R n_m for some n and m. Our goal is to show that the averaging 
function avr maps H± to 0. 

Let be an element of H±. This means that * r ) = for all G H T . 
The inner product (, } is invariant under the action of SO(A), so for every g 6 L, 

o vfr .?) = 0. But "fr is, by definition, invariant under T, so ° 9 = ^r, 
hence ° g, ^r) = for all *r 6 Hy- In other words, $iopg for all g 6 L. 
This implies that avr^jj G But avr(^j_) is also an element of Hy- Because 

Hy n H± = 0, this implies avr^jj = 0, as required. QED 

3.3.3. Interpreting the norm The squared norm |avr(^E'/3)| 2 provides a statistical 
estimate of the dimension of the L-invariant eigenspace. 

Proposition II. For each (3, the expected value of the squared norm (avr(^ r l a), avr('I , / 3)) 
is the dimension of the T -invariant eigenspace Hy- 



Proof. The comments in Step 1 of the algorithm in Section 3.3.1 show that the 
expected squared norm of the original i&p (before averaging) tells the dimension of 
the full eigenspace: 



2 



E{{^ f} ^p))=E{^$ >im ) = (3 

Choose an orthonormal basis {3^} for the eigenspace such that the basis vectors 
{3^1, ■ ■ ■ ,3m} span the T-invariant eigenspace Hy while the remaining basis vectors 
{3m+i, ■ ■ ■ , 3n} he orthogonal to Hy- Because the basis {3^} is orthonormal, the 
distribution exp(— r 2 /2) factors as a product 

n 

cxp(-r 2 /2)=n cx P(-C 2 /2)- 

i=l 

Proposition I implies that, relative to this basis, the averaging operator avp() preserves 
the first m coordinates while collapsing the remaining n — m coordinates to zero. Thus 
the distribution of avr(^/3) is given by the restricted product 

m 

nexp(-C 2 /2) = exp(-r 2 /2) 
i=i 

where r 2 = J^j=i d i an d the above reasoning now implies that 

J E((avr(* / 3),avr(* / 3)}) = m, 

as required. QED 

Proposition II remains valid even if we do not explicitly construct the basis {3^}- 
We may instead compute the squared norm (avr( v I / /3), avr( v I / /3)) by sampling points. 
Taking 32 random eigenmodes *S?p and for each one evaluating (avr(*/3), avri^p)) 
at 32 random points yields the following rough estimates for the dimension of the 
eigenspace (PDS designates the Poincare Dodecahedral Space of order 120): 
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k 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


S 3 


9.1 


17.1 


22.2 


40.0 


44.1 


60.7 


83.9 


106.4 


129.5 


144.6 


163.9 


p3 


8.5 


0.0 


23.9 


0.0 


47.2 


0.0 


80.5 


0.0 


120.2 


0.0 


165.8 


L(5,2) 


1.4 


3.6 


3.9 


6.6 


7.9 


12.7 


16.5 


19.0 


26.9 


25.1 


31.9 


PDS 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


11.8 



The above computations took only a few minutes on a 300 MHz desktop computer. 
Increasing the number of random eigenmodes and the number of sampling points 
would increase the accuracy of the results at the expense of a longer computation 
time. 



3.3.4- Computational complexity The projection method is reasonably fast. The 
choice of the Cf3im in Step 1 requires only O{0 2 ) time, and is completed almost 
instantaneously on a desktop PC. Thereafter the evaluation of avr(*/3)(a;) for each 
point i£S 3 takes |T| times as long as evaluating the underlying random eigenmode 
Vp(x) of S 3 at the same point x. In other words, using the projection method, an 
eigenmode of S 3 /T is |T| times as expensive to compute as an eigenmode of S 3 . More 
precisely, the time to evaluate ^ p{x) grows as |T| /3 3 , because for a fixed value of (3 
the eigenspace has dimension /3 2 , meaning that there are 1 terms to evaluate, each 
of which requires 0{(3) steps. In practice we evaluated the eigenmodes of S 3 using 
the toroidal coordinates method of p2|, but in principle the same runtime could be 
obtained using spherical coordinates. 

4. Analytical solutions for lens and prism spaces 

Besides the numerical methods presented above, there are special cases for which the 
eigenfunctions can obtained analytically [H. The results are indeed of importance 
to test the accuracy of our numerical computations. The method is based on the use 
of torus coordinates and applies to lens L(p,q) and prism S 3 / D* m spaces. We briefly 
recall the main points and results of [H . 

4-1- Summary of the general method 

The key of the method is to choose a coordinate system that respects the holonomy 
group r. We introduce the coordinates in R 4 , (x, y, z, w) by 

x = cosx' cos 8' 
y = cosx' sin 9' 
z = sinx' cos(p' 

w = sinx' shop' (23) 

so that the equation of the 3-sphere is simply x 2 + y 2 + z 2 + w 2 — 1 . Note that they 
are different from the 4-dimensional coordinates (^) introduced previously and that 
now the intrinsic coordinates have to range as 

< x' < < 6' < 2tt < <p' < 2tt. (24) 

For each fixed value of x' G [0, 7r/2] , the 9' and ip' coordinates sweep out a torus. Taken 
together, these tori almost fill S 3 . The exceptions occur at the endpoints x' = and 
x' = 7r/2, where the stack of tori collapses to the circles x 2 + y 2 = 1 and z 2 + w 2 = 1, 
respectively. 
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Identifying the eigenmodes of S 3 /T and the T-invariant eigenmodes of S 3 , as 
explained in our introductory remark, the eigenmodes of a lens or prism space are 
given by Z p -invariant or -D^-invariant eigenmodes of S 3 . 

An elementary construction ]42| ] shows that for each wavenumber k, with 
eigenvalue k(k + 2), the corresponding eigenspace of S 3 is spanned by the basis 

B k = { Qu m | \£\ + \m\ < k and \£\ + \m\=k (mod 2) } (25) 

where 

Qkim = cosW X ' sinl™! X ' Pi m| ' W (cos2 X ') 



P. 



\m\,\t\ 



x (cos \£\6' o 
being the Jacobi polynomial 



sin \£\0 ) x (cos |m|(^>' or sin|m|</?') (26) 



r'r"" W = i i ( H , + ") (^'J + _f) (x + D' (x - 1)"- 



and cos (resp. sin \£\9') being used when i > (resp. £ < 0), and similarly for the 



choice of cos or sin \ m\tp' . 



It is straightforward to see how the generating isometry of a lens space L(p,q), 
given in rectangular coordinates by 



/ cos27r/p — sin27r/p 

sin 2n/p cos 2ir /p 

\ 

or in toroidal coordinates by 






cos 2nq/p 
sin 2nq/p 



\ 



— sin 2nq/p 
cos2nq/p J 



(27) 



Y 



X 
9' 
i 



2-k/p 
- 2nq/p, 



(28) 



acts on the Qkim- The eigenmodes of L(p,q) comprise the fixed point set of this 
action. A set of simple numerical conditions tells how to select an orthogonal basis for 
this fixed point set, essentially as a subset of the basis Specifically, the eigenbasis 
for L(p,q) includes 



k(' in 



— < 



Qfcoo 
Qk±eo 

Qk0±m 

Qktm + Qk 



-Qke 



Qki 



V2 



„ + Q fc 



V2 



V2 



always 

ffi£ = 0(modp) 
iff qm = 0(modp) 
iff £ = qm (raodp) 
iff £ = — qm(modp) 



(29) 



For details as well as for the explicit form of the eigenmodes of prism spaces, please 
see Ref. [ ff2) . A similar analysis yields an explicit eigenbasis for a prism space. The 
Qkim are already mutually orthogonal, so after normalizing them to unit length we 
may use the above basis to construct unbiased random eigenmodes of a lens or prism 
space, with wavenumber k. 
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4-2. Extracting the coefficients S,p a im 

The previous analysis gives the decomposition of the eigenmodes on the basis Qpt Tl 



as 



*/3,s = ^2r}p,stmQfUm(x',0',<P') (30) 

but what we need are the coefficients £,p. s ern of the decomposition on the basis ypem 
as given in Eq. (|TT1). 

ypim and Qpi m are two basis of dimension (k + l) 2 for each [3 = k + 1, so up to 
a change of coordinates between toroidal and spherical coordinates, we can write 

Qf3lm = ^ Oipu'mm* ypl'm'- (31) 

Note that (£,m) and (£',m') do not vary in the same range since < I' < j3 — 1, 
\m'\ < t' and \t\ + \m\ < f3, \£\ + \m\ = f3 — 1 (mod2). The coefficients apu' mm ' are 
explicitly given by 

Q{3££' mm' / Qi3tm{x i 0', tf'Wpe'm'iXi 8, <f) sin 2 xsin9dxd6dip (32) 

where (x' , 0' , if') are functions of (x, 0. if ). It can be checked from (|[) and (|23| ) that 
x' = arccos [cos x cos 9] (33) 
cos x sin 9 ^ 



yl — cos 2 x cos 2 6* 

^ = V . (35) 

It can be checked that the two basis differ by more than just a change of coordinates. 
It follows that for the particular case of lens and prism spaces, the computation goes 
as follows: 

(i) Determine the coefficients r\p, s im as given in Ref. Wh \ thus mainly a tedious 
book-keeping operation. 

(ii) Compute the coefficients of the change of basis ([[l]); this has to be performed 
once for all for all spaces. 

(iii) The required coefficients are obtained by a matrix multiplication 

(,fi,stm = ^ Clpit'mm'np^l'm' ■ (36) 



4-3. The simplest example 

For the projective space, S 3 /Z2, there is only one generator which brings any point 
to its antipodal point. It follows that the eigenfunctions are just the average of the 
spherical harmonics evaluated with their antipodal equivalent. This can be sorted out 
analytically and one easily obtains that 

C = \ [1 - i-lf] S. t(tm) (37) 

where the function <5 s .(£ m ) = 1 if the label s can be identified with (£, m) and is zero 
otherwise. We recover the results from figure |l|, that is that there is no modes for (3 
even and that the dimension of the eigenspace for [3 odd is equal to the one of S 3 , 
that is (k + l) 2 . 
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5. Numerical results 

An interesting property is the distribution of the number of modes per wavenumber 
interval. In hyperbolic manifolds, the number N(< q) of modes smaller than q is well 
described by the Weyl asymptotic formula 



Vo\(H 3 /T) 
6tt 2 



iV(<g)~ --TV'V -l) 372 (38) 



for q 3> 1. In the case of the 3-sphere, q 2 = k(k+2) with multiplicity mult(g) = (fc+1) 2 
so that 

N ^ {<q) = ± (p+ lf = (fc + l)(fc + 2)(2 fc + 3)-30 

6 

and the Weyl formula, which applies also to the spherical manifolds, tells us that 

iv [rl ~7v [s3l /|r|. (40) 

We computed the eigenmodes and eigenfunctions of some spherical spaces with 
the different methods described above. First, it was checked that the spectra agree 
with the theoretical ones in the cases described in table [l]. In the particular case of 
lens and prism spaces, the eigenmodes and eigenfunctions agree with the ones obtained 
analytically in Ref. [E2|. 

The computational time of the averaging method is experimentally found to be 
proportional to |r|/? 4 68 , the typical running time being less than 10 seconds on a 
desktop computer for (3 < 15. 

In figures || and 5, we present the two examples of the lowest modes of S 3 /D\ 
and S 3 /Z%. In figure 6 we depict one of the ten eigenmodes of S 3 / D% with k = 4 for 
different values of the radial coordinate %. 

6. Some cosmological implications 

The goal of this section is not to compute the CMB anisotropics in details (this task 
will be delt with in a follow-up article jl3) ) , but to give estimate of the expected effects 
on large angular scales. 

The evolution of the scale factor, a, of the universe is dictated by the Friedmann 
equation that can be recast under the form 

TL \ 2 

— = n roX - 2 + n^ar 1 + n Ao x 2 + (1 - n ro - n mo - OO (4i) 

v lUlJ 

where H. = a' /a, a prime denoting a derivative with respect to the conformal time 77. 
We have introduce x = a/ag = (1 + z)~ , z being the redshift and ao the value of the 
scale factor today. The density parameters are defined by 

_ np ro _ np mo _ np\ g 

r ° — m — o ff 2 ' " A o — ott2 

O-Hq "J^o 

where n = 8ttG and Hq = Tio/aQ. In terms of these quantities, the physical curvature 
radius today is given by 

i?P hys ee a R Co = 4 1 (43) 

We can choose ao to be the physical curvature radius today, i.e. ao = Raj*' w ^ c ^ 
amounts to choosing the units on the comoving sphere such that Rc = 1, hence 
determining the value of the constant ao = c Oa + ^m — 1 /Hq. 
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Figure 4. The ten eigenmodes of S /-DJ f° r & = 4 in Hammcr-Aitoff projection 
for x = 7r/2. Note that some are identical but only on the equator and it can be 
checked that these modes are indeed linearly independent. 




6.1. Generalities 

When studying the CMB anisotropics, one has to go beyond the homogeneous and 
isotropic description of our universe and need to condider a perturbed spacetime with 
metric 

ds 2 = a 2 (77) [-(1 + 2<S>)drf + (1 - 2%)i ij Ax i dx j ] (44) 

where we consider only scalar modes and working in longitudinal gauge. Vector and 
tensor modes would have to be added for a complete description (including for instance 
gravitational waves, but are negligible on large angular scales. On these scales, we 
can neglect the effect of the anisotropic pressure so that *f> — <£> and the effect of the 
radiation between the last scattering surface and today. 

Under these assumptions, the temperature fluctuation in a direction n can be 
related to the gravitational potential <£> by, again in the particular case of 

adiabatic initial perturbations, 

— (n) = -$[r? LSS , (r? - ?7 LS s)n] + 2 jf 1 ^ " U r, (45) 

where ?y LSS and 770 are the value of the conformal cosmic time at the emission of 
the photon (last scattering surface) and at the reception (observer). According to the 
standard nomenclature, we will refer to the first term as the ordinary Sachs- Wolfe term 
(OSW) and to the second as the integrated Sachs- Wolfe term (ISW). The temperature 
angular correlation function is then defined by 

C(9) = /S ni )S,i 2 )\ . (46) 

\ / cos 9— ni .ri2 

Decomposing the temperature fluctuation on the spherical harmonics as 

f(n)=E E a tm Y t m{6,<t>) (47) 



Eigenmodes of 3- dimensional spherical spaces and their application to cosmology 19 




Figure 6. A mode of S3/D^ with k = 4 for 10 different values of x = iir/20 
(i = 1..10). It must be read left to right, top to bottom. 
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the coefficients Cg of the development of C(0) on Legendre polynomials are given by 

i 

{2i + l)C t = £ (<**A>- (48) 
m=—t 

To compute these quantities, one needs to determine the gravitational potential 
$. Its evolution is dictated J22I by the equation 

3>" + 3H(l + c 2 s )& - c 2 A 2 <I> + [2H' + (1 + 3c 2 )(H 2 - K)]G> = 0. (49) 

This relation strictly holds only for initial adiabatic perturbations as predicted by 
most of the inflationary scenarios, c 2 = P' j p' is the sound speed and is given by 
c s = Pr/(Pm + 4/3p r )/3. After decomposing $ on the eigenmodes as 

$(ry,x) = ^$ /3s ^( x ), (50) 

one can easily show that if the universe is matter dominated between the last scattering 
epoch and today then 

= F(r?)$ (q) (51) 

where we use the notation that q = (/3, s). $o is the value of the gravitational potential 
at the beginning of the matter era, but since in the early universe the curvature 
term is negligible and the dynamics is dominated by the radiation (so that a oc rf) 
it can be shown that, for long wavelengths (krj < 1), the non decaying mode of the 
gravitational potential is constant so that $o is in fact the primordial gravitational 
potential. Inflationary theories predict that it is a Gaussian field, and that all modes 
are independent, with power spectrum |fl6f 

2tt 2 

($ (q)*5(q')> = W 2 _ K f ^W - P')Sss'- (52) 

where K = 0, ±1 is the sign of the curvature. In the case of a scale invariant Harrison- 
Zel'dovich spectrum, V<s> oc [3°. 

Now, inserting the decomposition ( |50|) with the solution (^lj) in Eq. ( |45| ) and 
decomposing the eigenmodes as in Eq. one finally gets that the coefficients of 

the development ( [47| ) are given by 

aim = $o(/3, s)^f3 tS i m Gpi (53) 

0,s 

where Gpi is defined by 

1 fLSS 

GjH = ^F(n, ss )R pe ( m - r? LSS ) + 2 / F'{n)R pt { m - n)dn. (54) 



'/n 



The correlation 



2ir 2 

(atm,at'm>} = X! Wfp _ K ^ ^P) G PZ G f3t'£,P-Mrn£,*f3,sl>mi ( 55 ) 

has non-zero off-diagonal terms which reflects the fact that there is a global anisotropy 
due to the non-trivial topology. In a simply-connected homogeneous and isotropic 
universe {atmO-i'm') = Cg8u>5 mm i . These off-diagonal terms are characteristic of the 
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non-Gaussianity induced by the topology. From the expression (55), we can extract 
the Ci which characterise only the isotropic part of the temperature distribution, as 

p,s,m ™^ ' 

Indeed, in the Euclidean and hyperbolic cases, the sum over j3 has to be replaced by 
an integral; in the spherical case, (3 > max(3,i? + 1). This result has to be compared 
to its covering space analog, obtained by setting ^p^sim — 1 when the index s can be 
assigned the value (£m) and zero otherwise, so that the sum over m gives (2£ + 1), 

so that there is an average effect with the ponderation |£/3 iS f m | 2 /(2£ + 1). 

The effect of the topology on large scales can thus be investigated, in a first step, 
by considering the index 

Ti n = |r|^. (58) 

At large angular scales, will oscillate due to the missing eigenmodes while it 
will converge to unity, mainly because of Eq. (|40|), on small angular scales where 
the topology becomes irrelevant. Indeed, this will allow to put constraints on some 
topologies but an unambiguous detection will have to use a full-sky CMB map. 

6.2. Example of the three-torus 

As a first example, let us consider a cubic 3-torus of comoving size L. In the simplest 
case in which J7a = 0, F is constant and 

^'Ew*-^' (59) 

k 

with k = 27r(m, Ti2, ns)/L, and for the universal covering space we have 

C f oc J ^jj{u) oc (60) 

known as the Sachs- Wolfe plateau. 

On this simple example, we can determine the cut-off £ cu t below which there is a 
suppression of the power spectrum. The first approach is to remind that the Bessel 
functions peak at k(rjo — rj LSS ) ~ £. Hence a mode k contributes maximally to the 
angle subtended by its corresponding scale at last scattering. In flat models, rjo is 
given @ by 

^-^=(i+mC 5 ). (61) 

It follows that the OSW term has a cut-off round 

gcut- - „ 4 ^ (l + lnCT- (62) 



r. 



"mo 

This is analogous to the estimate by Inoue 
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Another approach, first introduced in p9| , 
which the maximum comoving wavelength A m; 

„ a L 



50[ , is to compute the angle & cut under 
£ = L at last scattering. It is given by 

(63) 



(l + z)d A (*) 

where dA = dL/(l + z) 2 is the angular distance and where the luminosity distance is 
given by 

du 



d h (z) = (l + z)^f 



(64) 



(l + u)E(u)_ 

In the case where fi A = 0, d L (z) = 2a (l+z)(l-(l + zy 1 / 2 )/H so that 6 cut ~ H L/2. 
We have introduced E = TL/TLq. Now, the £ th Legendre polynomial is a polynomial 
of degree i in cos# and has t zeros in [—1, 1] (or 2£ zeros in [— tt,tt] if working in 9) 
with approximatively the same spacing. We can estimate that ~ 7r/£ and thus that 
the cut in the OSW contribution is expected to be round 



2?rc 



f 



du 



{l + u)E(u)_ 

the factor 2 arising from the fact that an oscillation corresponds to 2 zero, 
corresponds exactly to the previous estimate (p2[). 



(65) 
This 



6.3. Spherical spaces 

In the spherical case we expect the same kind of effects, i.e. a suppression of the large 
scale ISW effect due of the existence of a maximal wavelength. The OSW term will 
be approximatively the same as the one computed in the flat case since our universe, 
even if closed, is still very flat. 

Let us first remark a crucial difference between Euclidean manifolds and spherical 
manifolds. For the former, as we have seen in the previous section, the smallest 
multipole is directly related to the size L of the fundamental polyhedron so that one 
cannot consider too small universes. For spherical universes, the situation is a priori 
different. As seen on the example of single-action manifolds, the value of the first 
non-zero eigenvalue does not depend on the order of the group, at least for cyclic and 
binary dihedral groups. 

Let us take the example of lens spaces, the first nonzero eigenvalue is indeed the 
same for all S 3 /Z m (m > 1), namely k = 2 and eigenvalue 2(2+2) = 8. However, the 
multiplicity is 3 for homogeneous lens spaces L{p 1 1), but only 1 for nonhomogeneous 
lens spaces L(p,q), as shown in Table 1 of [^2|. This can be understood by the fact 
that the space is becoming smaller and smaller in only one direction. In perpendicular 
directions the space remains large. The waves do not have distinct peaks (like 
mountaintops found in nature) but rather have ridges (perfect horizontal ridges, which 
are never found in nature). First imagine a wave in S 3 , as follows: the set of maxima 
is a great circle, while the set of minima is a complementary great circle. As time 
passes, the wave goes up and down, so what is the top at a time t = becomes the 
bottom at time t = 1/2 say, and vice versa, returning to its original position at time 
t = 1. Midway between the "top ridge" and the "bottom ridge" is a torus which 



remains at height for all times t. In toroidal coordinates (23), one ridge is the circle 
at x — 0, the other ridge is the circle at \' = n/2, and fixed torus lies at x' — 7r /4- 
This wave is preserved by all corkscrew motions along the natural axes, so it projects 



Eigenmodes of 3- dimensional spherical spaces and their application to cosmology 23 



down to an eigenmode of all lens spaces L(p,q). In the notations of Section ^. it is 
given in toroidal coordinates as 

*20o(Y, 0', 0') = sin 2 X ' - cos 2 X ' = ~ cos2 X ' 

thus verifying that it is a wave as described above, and in particular with no 
dependence on 9' or <fi' . In the special case of a homogeneous lens space L(p, 1), 
we get two more eigenmodes 

sin(2x') cos(0' - cj>') and sin(2x') sin(0' - <t>'). 

These modes are constant along helices (where 9' — <j>' is constant), so they are 
eigenmodes of all homogeneous lens spaces, whose holonomies are Clifford translations, 
but not eigenmodes of nonhomogeneous lens spaces. 

We thus expect the effect of the cut-off in the Sachs- Wolfe plateau to be milder 
than for Euclidean spaces. But, on the other hand we expect to have a more irregular 
Sachs- Wolfe plateau due of the fact that some wavelengths are missing from the 
spectrum. Another observational consequences missed by the angular power spectrum 
is a global large scale anisotropy that, at least, is expected for non-homogeneous lens 
spaces. 

To illustrate this, consider the simplest example S 3 jZ^ for which half of the 
modes have disappeared so that k = 2p, p being an integer. It follows that for I ~ 2p, 

[^2 ] [ "-^ ^ ] [ ^2 ] 

2 ~ C t and for I ~ (2p + 1), C e 2 ~ 0. We thus expect that the Ci curve 

oscillates around C e /2 with a frequency of order I ~ rj . 

Using the coefficients ([37]) for the projective space, we plot in figure [?], the function 
Ti for I < 20 for different curvature radius. We do not include the cosmic variance. 
This confirms the previous semi-analytical analysis and is in agreement with the 

numerical computations performed for S 3 /Z 2 ]2lf . When \Ct — 1| <C 1, C^ 2 ' ~ 1 
because x LS s < ^l^i the topological scale. When x LSS ~ 7r /2 ) both antipodal points 
on the last scattering surface are close to the equator so that the angular correlation 
function in opposite directions is expected to be higher. The case in which x LSS > tt/2 
is even more intricate because the geodesies are warping around the universe more 
than once. 

Such features are general to all spherical spaces, but the higher the order of the 
group, the larger the minimal curvature radius to get a topological signal. Since the 
size of the manifold decreases with the order of the group, there will always exist 
potentially detectable topologies even for spaces very close to flatness. 

Note also that in the case of the projective space, 

which can be understood if one remembers that there is no breaking of global isotropy 
and homogeneity for the projective plane. This is the only exception for which one 
can have a non-trivial topology and no preferred direction. 

7. Conclusion 

In this article, we have investigated in details the structure of the eigenmodes of 
the Laplacian operator in spherical spaces. A series of analytical and mathematical 
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Figure 7. The T e for I < 20 for S 3 /Z 2 - From top to bottom and left to right the 
curvature increases as Q = 1.1, 1.25, 1.5, 1.7, 1.8, 1.9 with A = 0. This corresponds 
respectively to x LS s = 0.5939,0.8990,1.1945,1.3528,1.4173,1.4746. For higher 
curvature, the topologicaal scale has more and more importance. The growth of 
the amplitude of the oscillations at small multipole is dur to the growth of the Ct 
due to the integrated Sachs-Wolfe effect. 
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results have been either reviewed or obtained and we have introduced various efficient 
numerical methods to compute them. These methods were compared together and to 
some analytical results. 

We have also investigated some cosmological consequences of this work par- 
ticularly concerning the large angular scales of the cosmic microwave background 
anisotropies. The effect of the topology in CMB calculation has been described and 
these results are now being included to a full Boltzmann code to simulate CMB maps 
with a topological signal [Q. As an example, we considered the simplest of all cases, 
that is the projective space. 



Appendix A. Eigenmodes of the Laplacian of homogeneous and isotropic 
simply connected three-dimensional spaces 

This appendix follows the work by Abbott and Schaeffer |5]| and Harrison |52| . Its 
goal is to summarize the derivation and explicit forms of the scalar harmonic functions 
solutions of the Helmoltz equations ([!]) . 

We rewrite the Friedmann metric (0) as 



ds 2 = -dt 2 



a 2 (r) 



dn 2 ) 



(A.l) 



with tr(r) = 1 + Kr 2 /4. By splitting radial and angular variables, it can be shown 
that the eigenmodes decompose as 

y e m (r, 6, 0) = R pe (r) Y em (0, cp). (A.2) 

Yi m (9, </>) are the spherical harmonics, related to associated Legendre polynomials P™ 

by 

' (2£ +!)(£- m)\] 1/2 



Y, 



P?(cos6) e im<s > 



Att(1 + m)\ 
and satisfy the relation 

WM) = (-ir^UM). 

The radial eigenfunctions Rpi(r) are solutions of the radial harmonic equation 



dR 



a 3 (r) d 
r 2 dr \o~(r) dr 



pt 



2 2 ( v g(£+i) 

q -a (r) — 



Rpi = 0. 



(A.3) 
(A.4) 
(A.5) 



In flat space (K — 0) the radial eigenfunctions are simply the spherical Bessel 
functions Rpi = ji(flr). 

The cases K = ±1 can be treated simultaneously by using the variable £ defined 

by 

sinC^ 172 -^- (A.6) 
cr(r) 

Note that £ is related to the dimensionless comoving radial distance in units of the 
curvature radius by £ = K 1 / 2 x- In terms of these variables the radial equation (A.5) 
takes the form 



1 



sin 2 £ d£ 



3; staT 



dR 



d£ 



sin £ 



R 



■pi 



0. 



(A.7) 
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Introducing the function Hgt(£) = Rpi(£) sin 1 / 2 ^) allows to solve the radial equation 
in terms of associated Legendre functions Pj^(cos£). 

For K = +1 (i.e. (3 2 = q 2 + lQ the radial eigenfunctions are given by 



Rptix) <* - 



-1/2-1 

smV - 1/2+P 



-P 



(cos x) 



with (3 > max(3, £ + 1) being an integer^. 

For K = — 1 (i.e. /3 2 = q 2 — 1) the radial eigenfunctions are given by 



Rpi(x) « — 



1 



-P 



-1/2- 
-1/2+i, 



^(coshx) 



(A., 



(A.9) 



sinh x 

and (3 can now take on any positive real value since there are no periodic boundary 
conditions to satisfy. 

We use the normalisation condition 

r 2 c\rc\£l 

y} tm yp't'm -T- = S(P - /?') &u> Smm> , (A.IO) 

so that the properly normalized functions take the following form 



N, 



pi 



1/2 



P 



-1/2-1 



Rpiix) = ' 



sinh X y - -i/2+i/9 

(2(3 2 /n) 1/2 M ((3 X ) 
( Mf3i\ 1/2 



(coshx) K = — 1 



k Vsinx/ 



P 



-1/2-e 

-1/2+P 



(cos x) 



K = 
K = +1 



(A.ll) 



with the two coefficients 

e e 
N PI = \{{I3 2 +n 2 ) M [M =\{{(3 2 -n 2 ). 



(A.12) 



n=0 n=0 

The radial eigenfunctions ( |A.ll ) differ from those determined by Abbott and 
Schaeffer Q by an overall factor (2{3 2 /ir)^ 1 / 2 due to the fact that they used the 
normalisation 

y*pi m yp'vm' ~^?~r = ~^32 S (P ~ Su ' 6mm '- (A.13) 

To finish, in the case of spherical space, the harmonic functions can be expressed in 
terms of the 4-dimensional coordinates (|^) as 

1/2 / . , \ V2 

I ' ! —I" I Ili III I ' 

yMrn{x) 



(2e + i)(e- m y. 



4tt(£ + my. 



M, 



1 M 



p 



-1/2-1 
■1/2+P 



(xo)P"t 



■1-2 



Xi + ix2 



yj\-x 2 j y/x\+. 



(A.14) 



Those expressions are of little value for numerical computation. There are two 
routes to compute numerically the eigenmodes. First, and as explained in Abbott 
and Schaeffer one can use a recursive relation between Rpt, Rpj-\ and i?^_2- 
Another efficient method |53j makes use of a WKB approximation. These two methods 
are complementary. 

+ We recall that /3 2 = q 2 + K. 

* It is well known that homogeneous harmonic polynomials of degree k on R 4 restricted to S 3 are 
eigenmodes of the Laplacian with eigenvalues k(k + 2). It follows that f3 = k — 1 is necessarily an 
integer. 
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